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^ ! Abstract 

We show that both small mixing in the quark sector and large mixing in the lepton sector can 

O ■ 

■ be obtained from a simple assumption of universality of Yukawa couplings and the right-handed 

■ neutrino Majorana mass matrix in leading order. We discuss conditions under which bi-large mixing 
in the lepton sector is achieved with a minimal amount of fine-tuning requirements for possible 

o 

<N 



models. From knowledge of the solar and atmospheric mixing angles we determine the allowed 
values of sin #13. If embedded into grand unified theories, the third generation Yukawa coupling 



£f~^ I unification is a generic feature while masses of the first two generations of charged fermions depend 

o 

on small perturbations. In the neutrino sector, the heavier two neutrinos are model dependent, 



while the mass of the lightest neutrino in this approach does not depend on perturbations in the 
leading order. The right-handed neutrino mass scale can be identified with the GUT scale in 
which case the mass of the lightest neutrino is given as (mf^/Moi/T) sin 2 #23 sin 2 #12 in the limit 
^ ■ sin #13 ~ 0. Discussing symmetries we make a connection with hierarchical models and show that 

the basis independent characteristic of this scenario is a strong dominance of the third generation 
right-handed neutrino, M 1 ,M 2 < 10" 4 M 3 , M 3 = M GUT . 
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I. INTRODUCTION 



The masses of three generations of fermions in the standard model are scattered in six 
orders of magnitude between the mass of the electron (0.5 MeV) and the mass of the top 
quark (175 GeV). Neutrino experiments suggest that masses of neutrinos (< 0.05 eV) are 
another seven orders of magnitude lighter than the electron. Furthermore, the mixing angles 
in the quark sector given by the Cabibbo-Kobayashi-Maskawa matrix Vqkm are small, while 
the mixing in the lepton sector is large. The solar neutrino data point to large mixing and 
the atmospheric neutrino data require close to maximal mixing. This is one of the most 
challenging puzzles in elementary particle physics. 

A lot of effort was made in order to understand the origin of hierarchy and mixing in 
fermion masses 1\ . Some of the most promising models are based on grand unified theories 
[GUTs] in which all fermions originate from just a few multiplets of GUT gauge symmetry 
group. Putting all particles of one generation into the same multiplet (as it is, for example, 
in SO (10)) offers a very attractive possibility that their Yukawa couplings are all equal at 
the GUT scale in a similar way to the well-established gauge coupling unification. This 
works quite well for the third generation of fermions but fails for the first two generations. 
Explaining the spectrum of the first two generations requires very specific assumptions about 
the underlying theory at the GUT scale. 

In the hierarchical approach, the unification of the third generation Yukawa couplings is 
typically the starting point and in the leading order only this coupling is generated. This is 
usually achieved by imposing family symmetries. Yukawa couplings of the first two families 
are generated in the process of family symmetry breaking. Although these models can 
reproduce the observed spectrum and mixing, the level of model building involved casts a 
shadow on third generation Yukawa unification itself. Why should the third generation be 
special when it is so easy to build models which do not unify the first two? Small mixing 
angles in this approach can be understood as a consequence of the large hierarchy in mass 
matrices. However, in order to obtain large mixing in the lepton sector it is often assumed 
that neutrinos are very different; the neutrino Yukawa matrix is not hierarchical or the right- 
handed Majorana mass matrix has a special form 

ELQ 

, or the neutrino sector is completely 

random Q]. 

Hierarchy in quark masses can also be understood within a democratic approach |4(, in 
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which all Yukawa couplings are identical in the leading order and all differences result from 
small departures from universality. However, in order to generate large mixing in the lepton 
sector, neutrinos are again assumed to be special. This time it is required that the neutrino 
mass matrix is not democratic (it is diagonal or it has another special form). This is either 
a by hand assumption p or it can be achieved in specific models p . 

In what follows we show that both small mixing in the quark sector and large mixing 
in the lepton sector can be obtained from a simple assumption of universality of Yukawa 
couplings. After introducing notation in Sec. m we start with a discussion of two families 
only, Sec. lHI[ since the mechanism which generates large mixing in the lepton sector is easier 
to follow in this case. Within a democratic approach, assuming the same universal form of 
all Yukawa matrices and the right-handed neutrino Majorana mass matrix (in the leading 
order), we identify a condition under which large mixing in the lepton sector is achieved. 
We show that the universal part of the resulting left-handed neutrino mass matrix is washed 
out due to the seesaw mechanism and the dominant contribution to neutrino masses comes 
from small departures from universality. This is what distinguishes quarks from leptons. 
The case of three families is discussed in detail in Sec. IIVI 

The virtue of this approach is that all mass matrices are treated in the same way, provid- 
ing a simple framework which can be easily embedded into more fundamental theories. If 
embedded into GUTs, the third generation Yukawa coupling unification is a generic feature 
in this approach, while the spectrum of the first two generations of quarks and charged lep- 
tons crucially depends on small perturbations. This becomes even more obvious in Sec. \V\ 
where symmetries of this framework are discussed and the connection with hierarchical mod- 
els is made. It is shown that the distinguishing (basis independent) feature of this approach 
is the dominance of the third generation right-handed neutrino mass. 

Masses of light neutrinos do not follow the same pattern as masses of quarks and charged 
leptons. The two heavier neutrinos are given in terms of perturbations and so are highly 
model dependent, while the mass of the lightest neutrino does not depend on details of a 
model in the leading order. If the right-handed neutrino scale is identified with the GUT 
scale, the mass of the lightest neutrino is predicted and it is related to the elements of the 
lepton mixing matrix. Avoiding the necessity of introducing an intermediate scale for right- 
handed neutrinos makes this framework very predictive. Finally, we conclude in Sec. IVII 

After the first version of this paper was finished []| it was brought to our attention that 
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a democratic form of both Yukawa matrices and the right-handed neutrino Majorana mass 
matrix was previously suggested in |8j]. Some of our results, namely the conditions under 
which large mixing in the lepton sector is achieved, coincide with the findings in Ref. |8f. 

tin 

For related studies, see also |9l llOlj. 
II. NOTATION 

The masses of quarks and leptons originate from Yukawa couplings of matter fields to 
one or more Higgs bosons. When Higgs fields acquire vacuum expectation values [vevs] , the 
Lagrangian containing mass terms of quarks and leptons can be written as: 

£ m = -Vff IA {Yf)i j f R j + h.c, f = u,d,e,u, (1) 

where vj result from vevs of Higgs fields which couple to the corresponding quark or lepton, 
fh Wr) represent left-handed (right-handed) fields, and Yf are Yukawa couplings, in general 
arbitrary 3x3 complex matrices in generation space represented by subscripts The 
Yukawa matrices can be diagonalized by bi-unitary transformations: 

Yf = UfY f Vj, (2) 

where Yf are diagonal matrices containing mass eigenvalues and Uf and Vf are unitary ma- 
trices. The mismatch between diagonalization of up-quark and down-quark Yukawa matrices 
appears in the charged current Lagrangian in the form of the CKM matrix: 

Vckm = U U U\. (3) 



The smallness of neutrino masses can be naturally explained by the seesaw mechanism |ll| 
which assumes the existence of Majorana masses for right-handed neutrinos: 

C UR = --v\M VR v R + h.c. (4) 

where Mr is a matrix in generation space. When right-handed neutrinos are integrated out 
we obtain a Majorana mass matrix for left-handed neutrinos: 

M Vh = -vXM-X T , (5) 
which can be diagonalized by a single unitary matrix: 

M VL = U VL M VL UT. (6) 



And finally, the lepton mixing matrix which appears in the charged current Lagrangian is 
given as: 

U = U e Ut L . (7) 



III. DEMOCRATIC MATRICES: TWO FAMILIES 

Let us start with only two generations and let us assume their Yukawa couplings are 
universal in the leading order: 



If the Yukawa matrices are exactly equal to XA//2, mass eigenvalues are {0,A/} and the 
diagonalization matrix is: 



Therefore one generation is massless in the leading order and the CKM matrix is the identity 
matrix as a consequence of the unitarity constraint Vckm = U U U\ = UjUj = 1. This is 
quite a good approximation to reality, taking into account that the only assumption so far 
is the universality of Yukawa couplings. 

Now let us parametrize the departure from universality by matrices Sf so that [21^ : 



Taking, for example, e/n = and e/12 = e/22 = e/ we find m s /mb — ed/4, m c /m t ~ e u /4 and 
V c b — (e<i — e u )/4 ~ m s /mb. The predicted value of \V c b\ in this case is 0.02 which is not so far 
from the desired value |V^| = 0.036. Clearly, relaxing the condition between elements of S, 
there is enough freedom to fit precisely both quark masses and CKM elements. For specific 
models see Ref. The generic consequence in this approach is that the value of |T4&| is 
proportional to the generated hierarchy in quark masses. The diagonalization matrices differ 
only by 0(e) from Uj. Therefore the resulting CKM matrix differs from the identity matrix 
by 0(e) although the mixing in both up and down sectors is close to maximal. 

Let us turn our attention to the lepton sector. If the only source of neutrino mass was 
the neutrino Yukawa matrix, the lepton mixing matrix would naturally be close to the 








(10) 
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identity matrix (in the same way as the CKM matrix), although it is reasonable to expect 
off-diagonal elements to be much larger than the corresponding CKM elements due to the 
fact that the hierarchy in both charged lepton and neutrino sectors is much smaller than 
in the quark sector (m U2 /m u:j > 0.16). However, close to maximal mixing would require 
fine-tuning [t v \2 ~ 1). Fortunately, we will see that the lepton mixing matrix can naturally 
be very different when considering right-handed neutrino Majorana mass and the seesaw 
mechanism. 

Let us make a similar assumption about the form of M VR as we made for Dirac Yukawa 
matrices: 



2 

It is useful to write the inverse of this matrix as: 



M„ R = l(I-n)M , K= [ rU ^ | . (11) 

Tl2 r 2 2 



M;' = ^(^+n), (12) 



M eff 



where 



and the effective right-handed neutrino mass scale is 





f-i 




|. *-| 


' r 22 


-ri 2 | 










{ -ri2 


ni J 



(13) 



M eff = ^(r-det^)Mo, (14) 

with 

2 

r = ^ jjjj = rn + r 22 - 2r 12 . (15) 

In the case when are much smaller compared to e U ij from Eqs. (jSJ) and (|12j) we get the 
left-handed neutrino Majorana mass matrix in the form: 



M„ - ™ 



UL AM, 



eff 



e' 2 e' e 
e'e e 2 



+ rX + Ofr^e 



(16) 



where ej, = e„n — e vV2 and = — 6^22- The special form of X causes that (X — £ V )X(X — 
8 V ) T = S U XS^ which corresponds to the first matrix in Eq. (|16jl. The second term in Eq. (fTKj) 
comes from X7ZX and the last term includes —EJKX — XTZS^ + S U TZS^ '. The eigenvalues of 
the first matrix in Eq. ()16|) are {0, e'J + e 2 ,} and the diagonalization matrix is given as: 

X 1 



U UL = VT+^ x = ^. (17) 



vi+x 2 \Zi+^ 2 
6 



In the limit x — > ±00 or x — > this matrix is either diagonal or off-diagonal. Assuming 
further that the charged lepton diagonalization matrix has approximately the form of Uj 
in Eq. we find that the lepton mixing matrix J7J) is approximately equal to t/j up to 
sign changes in different elements corresponding to different limits. These situations lead to 
maximal mixing in the lepton sector and they occur when one of en and €22 either dominates 
or is close to £12. Both situations are reasonable to assume. No particular fine tuning is 
necessary. Close to maximal mixing, U\2 ~ l/v2± 0.1, will be achieved every time when 
\x\ > 7 or \x\ < 1/7. 

The second term in Eq. (fi7)|) lifts the first eigenvalue. In situations when large mixing in 
the lepton sector is generated, the masses of left-handed neutrinos are given as A 2 u 2 /(2rMo) x 
{r, e 2 }, where e 2 = max(e 2 , e'^) and we assume r < e 2 and r ^> r 4 2 - (in other words there 
is no cancellation between r^- which would make r smaller than higher order terms 22j). 
An interesting consequence of this approach is a very robust prediction for the mass of the 
lightest neutrino: 

m vx = A> 2 /(2M ), (18) 

which is given by the universal Yukawa coupling \ v and the overall right-handed neutrino 
mass scale M . It does not depend either on details of the Yukawa matrix (e vi j) or details 
of the right-handed neutrino Majorana mass matrix (r^). We discuss this more in the case 
of three families. 

Before we discuss three families, let us summarize why it was possible to obtain large 
mixing in the lepton sector. The form of X given in Eq. (fT3|) plays the crucial role. The 
sum of elements of this matrix in every row and column is zero. As a result, XX = and so 
the Is from the neutrino Yukawa matrix wash out, leaving products of e^j as the dominant 
contributions to the left-handed neutrino Majorana mass matrix after the seesaw. If the 
resulting M UL is hierarchical (and it is if we assume that S u is similar to the perturbation 
matrices for other fermions), the lepton mixing matrix will be dominated by U e . Perhaps 
the only nontrivial assumption is that r < e 2 . However, this assumption does not require a 
departure from a democratic approach. Quite the contrary, it just means that M UR has to 
be somewhat more democratic than Dirac Yukawa matrices. We will see that the matrix X 
has the same property in the case of three families although the form of this matrix and the 
resulting left-handed neutrino Majorana mass matrix is much more complicated. 
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IV. DEMOCRATIC MATRICES: THREE FAMILIES 



Let us now assume that all three generations are indistinguishable in leading order and 
so Yukawa couplings are given as: 



( e f f 

e /n e /i2 e/i3 



(19) 



e/12 e/22 e/23 
y 6/13 e/23 e/33 J 

where we use the same symbol X for the 3x3 matrix with unit elements as we did in 
the 2x2 case, and similarly we parametrize the departure from universality by matrices 
£f. If Yukawa matrices were equal to ZA//3, then mass eigenvalues are {0, 0, A/} and the 
diagonalization matrix is: 



(20) 



As a consequence of degenerate zero eigenvalues the first two rows of this matrix are not 
uniquely specified and are model dependent (S has to be taken into account). They can be 
replaced by any of their linear combinations and the corresponding orthogonal combination, 
which is accounted for by the first matrix which rotates the first two rows. As a result, the 
CKM matrix is not the identity matrix in the leading order as it was in the case of two 





^ cos 9j 


sin 6x 


o\ 


/ i 

V2 


i 

V2 


o\ 




— sin 9 j 


COS 6*x 





1 


1 

V6 


2 




{ 
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V 7s 


1 


73 / 



23|. 



families, but rather a unitary matrix with an arbitrary 1-2 element 

The Majorana mass matrix for right-handed neutrinos is parametrized in a similar way 
as before: 



M, 



(J -ft) M , ft 



^ T\\ 7*12 ri3 ^ 

r\2 r 2 2 r 23 
ri3 r 23 r 33 



(21) 



The inverse of this matrix is given as: 



V R 



M, 



eff 



J + ft 



(22) 



where X can be written as 



1 



/ 



-n n - r 3 +r 3 
T\ - r 3 — n + 2r 3 - r 2 — r 3 + r 2 
+^3 -J"3 + r 2 -r 2 



(23) 
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with 



ri = r 22 + r 33 - 2r 23 , 
^2 = 7*u + r 22 - 2r 12 , 
r 3 = r 22 + r 13 - r 12 -r 23 . 



(24) 
(25) 
(26) 



From the form of X it is easy to see the correspondence with the 2x2 case, since it can also 
be written as: 



X = riXi + r 2 X 2 + r 3 X 3 + X. 



rT 

-3 ; > 



where 



Xi 



V 



1 -1 




Xo EE 



/ 




-1 1 
1 -1 



X 3 EE 



/ 



V 



0-1 1 \ 
1 -1 
Oy 



(27) 



(28) 



The matrix 1Z is proportional to the inverse of the matrix 1Z. Its elements are cofactors of 
the corresponding elements of TZ: 



And finally, M e tt is given as 



1 



M eff = -{r-detK)M , 

where 

3 

r = ^2 = rir2 _ r 3- 
Let us again assume that are much smaller than e V y. In this case we get: 



(29) 



(30) 



(31) 



M, 



9M 



M + rl + 0(n i3 e ul3 ) 



(32) 



where M. = E V XE J is a 3 x 3 equivalent to the first matrix in Eq. ()16j). As in the case of 
two families the matrix M. depends on differences between e U ij. It is useful to define the 
perturbation vectors: 



g — £v2 — &u3, 



(33) 
(34) 
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where £ V i is the i-th column of the perturbation matrix E v . The matrix Ai can be written 
as: 

M = -r l (e.e T ) - r 2 (g.g T ) - r 3 (e.g T + g.e T ). (35) 

From this expression it is easy to see that M. Vq = for v J_ e, Vq J_ g. Therefore, the 
eigenvector corresponding to the zero eigenvalue is given as: 

ex q . . 

v = p—^r. 36 
|e x g\ 

The heavy two eigenvalues can be written as: 



m ± = Ml±yi-^|. (37) 



with 



where 



t = p% + p 2 + 2p 3 cos a, (38) 
= (P1P2 - Pi) sin2 «> (39) 

Pi = — | e | 2 , p 2 = -r 2 |(?| 2 , p 3 = -r 3 | e | | g\, (40) 



and a is the angle between e and g. The eigenvectors corresponding to these two eigenstates 
are given as two orthogonal linear combinations of e and g: 

v± = a± e + b± g, (41) 

where a±, b± can be written in terms of pi, p 2 , p 3 , and cos a. Before we proceed further let 
us summarize the current status of neutrino masses and mixing. 



A. Experimental results 

A global analysis of neutrino oscillation data 
mass-squared differences: 



Q gi 



12, hj gives the best fit to the neutrino 



Ami, = m 2 U2 - m 2 Ul ~ 6.9 x 10" 5 eV 2 , (42) 
Am 2 atm = ml - ml ~ 2.6 x 10~ 3 eV 2 , (43) 



10 



and mixing angles: 

sin 2 # 12 = sin 2 ^ oi = 0.30, (44) 
sin 2 23 ee sin 2 ^ m = 0.52. (45) 

The 3o" ranges for mixing angles are: 

0.23 < sin 2 # soi <0.39, (46) 
0.31 < sin 2 # atm < 0.72, (47) 

and the 3a upper bound on the third mixing angle is: 

sin 2 #i3 < 0.054. (48) 

In the case m Ul <C m„ 2 , m U3 we can interpret these results as: 



m U2 ~ yAm 2 ^ ~ 8.3 x 10~ 3 eV, (49) 
m U3 ~ y/Am 2 atm ~ 5.1 x 10~ 2 eV, (50) 



and due to sin 2 13 ~ 0, the mixing angles are related to the elements of the lepton mixing 
matrix in the following simple way: 

sin 2 # 12 ~ \ U 12 \ 2 , (51) 
sin 2 # 23 ~ \ U 23 \ 2 . (52) 

B. Back to three families 

From the experimental results above we see that the lepton mixing matrix is charac- 
terized by a very small 1-3 mixing and close to maximal 2-3 mixing. The lepton mixing 
matrix originates from both the charged lepton diagonalization matrix U e and the neutrino 
diagonalization matrix U„ L . In a democratic approach the charged lepton diagonalization 
matrix already contains large mixing angles. If also U UL contains large mixing angles, it 
would require a conspiracy between elements of all perturbation matrices in order to achieve 
sin 2 # 13 ~ and sin 2 #23 — 0.5. Of course, such a conspiracy might naturally occur in some 
models. 

In order to avoid any exact relations between elements of £ e , £„, and 1Z the simplest way 
to proceed is to assume that the perturbation matrix £ v introduces the minimal amount of 
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FIG. 1: cos U as a function of z, see Eq. 1)54(1 . 

mixing into the lepton mixing matrix. This happens when perturbation vectors e and g are 
dominated by a single element ~ (0, 0, 1) T , (0, 1, 0) T , or (1, 0, 0) T with e J_ g. In this case 
U VL is given as: 

' 1 

cos 9 U — sin 6 V 
sm9 u cos9 v 

where P is a permutation matrix which interchanges rows of the identity matrix depending 
on the choice of perturbation vectors e and g. Note that, since cos a = 0, the eigenvalues 
and cos^ are given in terms of 3 parameters: pi, p 2 and p 3 . Namely, in Eq. (JBTj) t = pi + p 2 
and d = pip 2 — pi, and in Eq. (J53J), 

1 



\ 



\ 



I 





(l 





o\ 


p 





1 







1° 





V 



(53) 



cos9 h 



^2 + 2z 2 - 2zVl + z 2 



(up to an overall sign), where 



z = 



P 2-P1 
2p 3 



(54) 



(55) 



This function is plotted in Fig. ^ We see that cos 6?^ varies fast for small \z\ and is almost 
constant for \z\ ^> 0. Therefore, the least-fine-tuned situations correspond to \z\ 3> which 
happens for |p 3 | max{|pi|, \p2\}- Actually, no strong hierarchy is necessary, cos 9 U is very 
close to 1 for z as small as 3. Note that, if we assume | e | <C | g\, which is reasonable to 
assume about all perturbation matrices in order to generate hierarchy between generations, 
|p 3 | max{|px|, |p 2 |} is satisfied for a huge variety of possible entries r^. 
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FIG. 2: Contours of constant sin 2 #23 (blue) as a function of 9 e and 6 U . Dark blue contours represent 
the central value and the blue shaded area corresponds to the 3a allowed region. Overlayed is the 
3o" region of sin 2 #13 (green stripes). The green shaded area represents the overlap of 3a allowed 
regions of sin 2 #23 and sin 2 #13. 

From Eqs. (j7J), (J2Uj) and (JoHj) we see that the most general form of the lepton mixing 



matrix in this case can be written as 
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— sin 6 V 


cos 6 V 


I 



(56) 



where the permutation matrix from Eq. (|53|) was absorbed in the redefinition of cos# e 
since the permutation matrix switches columns of the second matrix above and this can be 
accounted for by rotation of the first two rows of that matrix. 

The lepton mixing matrix in this simplest scheme is given in terms of two parameters 9 e 
and 9 Vl and so it is not trivial that the resulting three mixing angles can be simultaneously 
within experimental bounds, although, from the suggestive form of the matrix in the middle 
of Eq. it might be guessed that it will happen for small 9 e and 9 V . In Fig. Elwe present 
contours of constant sin 2 #23 (blue) as a function of 9 e and 9 U . Dark blue contours represent 
the central value (Eq. (|15j)) and the blue shaded area corresponds to the 3a allowed region 
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71 

3tt/4 
O v tt/2 
tt/4 

tt/4 tt/2 3tt/4 7t 

FIG. 3: Contours of constant sin 2 9\2 (yellow) as a function of 9 e and B v . Brown contours represent 
the central value and the yellow shaded area corresponds to the 3a allowed region. Overlayed is 
the 3a region of sin 2 #13 (green stripes). The pink shaded area represents the overlap of 3a allowed 
regions of sin 2 #12 and sin 2 #13. 

(Eq. ([47)1 ). Green stripes represent the 3a allowed region of sin 2 # 13 given in Eq. (|48[) and the 
green area is the overlap of 3a allowed regions of sin 2 6*23 and sin 2 6*13. We see that it is not 
particularly difficult to achieve a very small #13 and close to maximal 823 ■ In the majority 
of cases, the value of sin 2 623 consistent with the limits on sin 2 613 is within or above the 3a 
allowed region. 

In Fig.Elwe present a similar plot for the solar mixing angle. Contours of constant sin 2 9 12 
are represented by yellow. Brown contours correspond to the central value (Eq. (|44jl ) and 
the yellow shaded area represents the 3a allowed region (Eq. ([46))) . Green stripes have the 
same meaning as in the Fig. Eland the overlap of 3a allowed regions of sin 2 9 12 and sin 2 #13 is 
given by the pink area. The value of sin 2 9 12 consistent with the limits on sin 2 # 13 is always 
within or above the 3a allowed region. 

However, the 3a allowed regions of all three mixing angles overlap only in two small 
regions as can be seen in Fig. HJ These regions are represented by magenta (note the 
periodicity of the picture; regions disconnected at the boundaries of the plot are not counted 




14 



71 

3tt/4 
O v tt/2 
tt/4 

tt/4 tt/2 3tt/4 7t 

FIG. 4: The 3cr allowed regions of sin 2 #23 (blue, green)), sin 2 9\2 (yellow, pink) and sin 2 #13 (green, 
pink) from Fig. and Fig. 01 over layed. The magenta regions represent their overlap. 

separately) . 

C. How small can sin 2 #13 be? 

Since we do not measure #13, it is interesting to ask what values of sin 2 # 13 can be achieved 
within this approach while satisfying the 3cr experimental bounds of sin 2 #23 and sin 2 # 12 . In 
Fig. 0] we see that the 3a regions of sin 2 #23 and sin 2 9 12 overlap in 4 places: two regions 
containing magenta areas (Region I), and the regions around points 9 e ~ 7r/2, 9 v ~ 57r/8 
and 9 e ~ 77r/8, 9 v ~ 7-7r/8 (Region II). Note that the two areas in the Region I predict sin #13 
with opposite signs, and the same applies to the Region II. The predicted values of sin 2 9 13 
from these regions are: 

0.008 < sin 2 #13 < 0.14 (Region I), (57) 
0.22 < sin 2 #i3 <0.66 (Region II). (58) 

Region I overlaps with the experimentally allowed region and it shows that the value of 
sin 2 #13 which can be accommodated in this approach can be as low as 0.008. Note that this 
minimal value of sin 2 # i3 corresponds to the maximal allowed values of sin 2 #23 and sin 2 # i2 . 
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On the other hand, the central values of sin 2 #23 and sin 2 9u correspond to sin 2 #13 near its 
present experimental upper bound. 



D. Consequences for models 



The magenta regions in Fig. 0] are not large. Nevertheless, the good news is that they are 
located around 9 V ~ 7r/2 or 9 V ~ 0, 7T, which corresponds to the least-fine-tuned possibilities 
discussed after Eq. (J55J). 

The only allowed values of 9 e are close to or 7r/4. The charged lepton diagonalization 
matrix corresponding to these two possibilities is close to the forms: 



for 9 e = 0, and 
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(59) 



V 



(60) 



2 2 

for 9 e = 7r/4 (up to overall sign changes in the rows). The first solution is obvious, since U e 
is very close to the observed lepton mixing matrix in this case. The second solution is not 
so obvious and its symmetric form is quite surprising. 

Perhaps the simplest models of this kind are those with a dominant 3-3 component in 
all perturbation matrices, £ ~ 1Z ~ diag(0, 0, e). This perturbation generates masses of the 
second family. In this case U e is very close to the one in Eq. (JST?|) and g ~ (0, 0, 1) T . There are 
many ways to introduce masses of the first family. Simple examples are £ ~ TZ ~ diag(0, 5, e) 
(this perturbation was suggested in 8]) or diag(5, 0, e), where 5 e (with different numerical 
values for different perturbation matrices). These situations have all the desired features: 
e _L g (approximately) and \z\ 3> 0. 

It is certainly remarkable that very simple forms of perturbation matrices, which on top of 
everything can be chosen to be the same for all mass matrices, lead to the observed pattern 
of fermion masses and mixing. This should come with a warning, however. The simple form 
of perturbation matrices does not guarantee that it is easy to obtain them naturally in some 
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models. Especially in models with family symmetries this is very complicated compared to 
the situation in hierarchical models. The reason is that in hierarchical models hierarchy is 
achieved by suppressing or forbidding some entries in mass matrices. This can be easily 
achieved by assigning different charges under a family symmetry to different particles. In 
democratic models however, family symmetries have to allow all entries in mass matrices 
and yet they have to account for small differences in specific elements of mass matrices. This 
is certainly not trivial to achieve. 



E. Mass of the lightest neutrino 

The mass of the lightest neutrino is lifted when the second term in Eq. is taken 
into account. Since we assume it is just a small correction to the first term, it can be 
treated as a perturbation. Adding this perturbation does not significantly affect the two 
heavy eigenvalues and the diagonalization matrix, but it is crucial for the lightest eigenvalue 
which is exactly zero in the limit when this term is ignored. In the case of non-degenerate 
eigenvalues the correction to eigenvalues rrii of a matrix Ai generated by a matrix SAi are 
given as: 

Srrii = uj 5A4ui, (61) 

where Ui are normalized eigenvectors. In our case SAi = rl (up to the overall factor 
in Eq. (JH2j)) and the eigenvector corresponding to the zero eigenvalue is vq ~ (1,0, 0) T . 
Therefore, 

m - = 9M^T (62) 
Since M e ff ~ rM /3, see Eq. (jSQj), the mass of the lightest neutrino again does not depend 

on details of a model in the leading order and is given as: 

m - = 3Mo ■ (63) 
This result is based on our assumption of the minimal amount of mixing coming from the 
neutrino diagonalization matrix. However, it is possible to make a prediction which does 
not depend on this assumption. 

Let us suppose that we do not know what the eigenvector corresponding to the lightest 
eigenvalue is. Due to the universal form of SAi, we have 

v^SMv = r£ 2 , (64) 
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where 

3 



1=1 



and so the mass of the lightest neutrino is given as: 

A 2 ?; 2 

m =^JL^\ (66) 
1 3M V ; 

In general £ can be anything between and v3- However, in order to satisfy bounds on 
lepton mixing angles £ cannot be arbitrary. The 3-1 element of the lepton mixing matrix is 
given by: 

and so 

"V = ff^, (68) 

Note that the 3rd row in U e is not model dependent unlike the first two rows are! It can 
receive only small corrections from the perturbation matrix. Finally, in the case of complex 
matrices, v T in Eq. ()64|) becomes v* and U 2 X in Eq. (jUHj) becomes | U T \ | 2 . 

Although we do not measure U T i, it is related to the observed mixing angles due to the 
unitarity of the lepton mixing matrix. In the case sin 9 13 ~ it is simply given by 

U T i = sin #23 sin #12. (69) 



A global analysis of neutrino oscillation data 



gives the 3a range: 



0.20 < I U T i I < 0.58. (70) 

The value of U T \ = 1/ Vs in which case Eq. (JfiHJl gives the same result as Eq. (jHSJ) is close to 
the upper limit. 

The masses of the two heavier neutrinos are given in terms of r^, e u ki, and so they are 
highly model dependent. Let us look at a simple example to get a feeling for typical values 
of perturbations which lead to the observed spectrum. Let us assume that the form of all 
perturbations is ~ diag(5, 0, e). In this case we have: | e | 2 = <5 2 , | g\ 2 = e 2 , r 1 = e r , r 2 = 5 r , 
and r3 = 0, from which we get: pi = e r 5 2 , p 2 = S r el, = 0, and r = r\T2- The neutrino 
masses are given by: 

AV A 2 t» 2 X 2 v 2 (5 2 e 2 



A 2 f 2 X 2 v 2 (S 2 e 2 l 

3M r 3M I d r e r J 



9M e// 
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Using Eq. (|55|) we get 

In order to have m U2 ,m U3 > m Ul we need <5 r < <5 2 and e r < e 2 , and so the hierarchy in the 
right-handed neutrino mass matrix has to be much larger than the hierarchy in the neutrino 
Yukawa matrix. This coincides with the assumption we had to make in order to achieve 
large lepton mixing, see Eq. (JH2J- 

In simple SO (10) type models X u v u = X v v u , in which case the lightest and the heaviest 
fermion of the standard model are connected through the relation in Eq. (|68|) where A 2 t> 2 
is replaced by m 2 op (actually, to be precise, X u = X u is a relation at the GUT scale and the 
effects of the renormalization group running between the GUT scale and the electroweak 
scale should be taken into account). This is a very pleasant feature since we can further 
identify M with the GUT scale, Mgut ~2x 10 16 GeV, in which case we get 

m vl = ^ | U Tl | 2 , (73) 

and predict the mass of the lightest neutrino to be between 5 x 10~ 5 eV and 5 x 1CT 4 eV 
depending on the value of U T \. 

From experimental values of m V2 and m U3 , given in Eqs. (J49)) and (J5U|) . and from the 
predicted mass of m ui , we see that e 2 /e r in Eq. (|72|) is of order 10 2 and <5 2 /5 r is an order 
of magnitude smaller. The largest e/ necessary to fit the masses of the second generation 
of charged fermions corresponds to e e ~ 0.28. If we take e„ ~ 0.1, we find e r ~ 10~ 4 and 
the mass of the second right-handed neutrino is Mi ~ (10~ 4 — 10 _5 )M3, where M3 = Mq = 
Mqut- This is a very rough estimate and it is not clear how to reasonably estimate M l5 
besides the relation Mi < M 2 . We conclude that the spectrum of the right-handed neutrinos 
consistent with bi-large mixing in our setup is Mi < M 2 < 10 M3 and M3 ~ Mgut- 



V. A SYMMETRY FOR THIRD GENERATION DOMINANCE 

The democratic Yukawa matrices are well motivated by S3L x S^r family symmetry 
However, the right-handed neutrino Majorana mass matrix is then constrained by the S^r 
symmetry only and the democratic form is not unique. The most general form of the right- 
handed neutrino mass matrix is given as a linear combination of the democratic matrix and 
the identity matrix. 
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It has been shown that the democratic forms of Yukawa matrices and the right-handed 
neutrino mass matrix are uniquely specified by imposing a Z% symmetry realized in the 
following way 



fu 
fm 



PijfRji 



(74) 
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where 
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uj 1 1 \ 

1 UJ 1 
1 1 UJ 

It is straightforward to check that this is indeed a Z 3 symmetry, and the proof that it is 
responsible for the democratic form of Yukawa matrices and the right-handed neutrino mass 
matrix can be found in Ref. [^]. Here we provide an alternative proof and a very simple 
understanding of this peculiar symmetry. 

Let us first note that a democratic matrix Mo (which represents both Yukawa matrices 
and the right-handed Majorana mass matrix): 

ti i A 

1 1 1 



M D = - 
3 



1 1 1 



(77) 



can be brought to a diagonal (hierarchical) form: 



^0 ^ 



M, 



H 




1 



(78) 



by a unitary transformation 



(79) 



M H = U X M D U%, 

where Uj is given in Eq. (|20j) . 

It is obvious that a Z 3 symmetry under which the first two generations of left-handed 
fermions have charge -1, the first two generations of right-handed fermions have charge +1 
(or, equivalently, charge conjugates of right-handed fermions have charge -1), and the third 
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generation of fermions has charge zero uniquely specifies the Mjj form of Yukawa matrices 
and the right-handed neutrino Majorana mass matrix. It follows from the fact that any 
coupling involving first and/or second generation fermions is simply forbidden. Therefore 
the only non-zero element of a mass matrix is the 3-3 element. This can be written in the 
form of transformations (J74j) and (|75|) with P replaced by Pjj' 



P, 



H 



* uj 0^ 
u 
1 



(80) 



where the subscript H indicates that it is a symmetry transformation specifying the hierar- 
chical form of mass matrices. Therefore we proved: 

Pjj M Pjj = M M = XMjj. (81) 

Now we can rotate this result to the democratic basis: 

U^PhUjU^MUjU^PhUj = U^MU 2 U%MU X = \U%M H U X , (82) 

and we obtain: 

PM'P — M' <=> M' — XMd, (83) 

where M' = U^MUj and 

P = U^PhUj. (84) 

Inserting Pjj from Eq. (|50jl and Uj from Eq. (|20|) it is straightforward to find that the form 
of the transformation matrix P which guarantees the democratic form of all mass matrices 
is exactly that of Eq. (J75j) . 

We see that the special form of P in Eq. (J76|) is just a simple Z 3 symmetry which allows 
only 3-3 elements in mass matrices rotated into the democratic basis. This shows that the 
approach in which the right-handed neutrino Majorana mass matrix and Yukawa matrices 
have in the leading order democratic form is equivalent to the hierarchical approach in which 
the 3-3 elements in Yukawa matrices and the right-handed neutrino mass matrix dominate. 
The results obtained in previous sections can be translated into results in corresponding hi- 
erarchical models Namely, the prediction for the mass of the lightest neutrino does not 
change, since mass eigenstates are not basis dependent. The necessary (and basis indepen- 
dent) requirement for achieving large mixing in this scheme is M l5 M 2 < 10 _4 M 3 , i.e., the 
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third generation right-handed neutrino has to dominate even more than the third generation 
of quarks and charged leptons. Thus third generation dominance is a suitable name for this 
scenario. 

Besides this Z 3 symmetry, the democratic mass matrices have also larger symmetries, 
like S3, for example, which can be used to specify perturbation matrices in the process of 
family symmetry breaking and thus distinguish between hierarchical and democratic starting 
point |26|. However, the complicated form of the Z3 symmetry in the democratic basis 
compared to the simple form in the hierarchical basis further amplifies the difficulties in 
constructing specific models as discussed at the end of Sec. IIVDI 

The fact that these mass matrices are motivated by a simple Z 3 symmetry is certainly 
a very pleasant feature. Unlike S^l x S3R, this Z 3 symmetry acts in the same way on all 
particles in each generation. Therefore this approach can be readily embedded into GUT 
models, like SO (10). 



VI. CONCLUSIONS 



We showed that both small mixing in the quark sector and large mixing in the lepton 
sector can be obtained from a simple assumption of universality of Yukawa couplings and the 
right-handed neutrino Majorana mass matrix in the leading order. We discussed conditions 
under which bi-large mixing in the lepton sector is achieved with a minimal amount of 
fine-tuning requirements for possible models. From knowledge of the solar and atmospheric 
mixing angles we determined the allowed values of sin # 13 . The central values of sin 2 #23 and 
sin 2 #12 predict sin 2 #13 near its present experimental upper bound, while it can be as small 
as 0.008 if both sin 2 #23 and sin 2 #12 are near their 3a upper bounds. 

We showed that this approach is equivalent to the hierarchical approach in which the 3-3 
elements in Yukawa matrices and the right-handed neutrino mass matrix dominate. The 
necessary (and basis independent) requirement for achieving large mixing in this scheme is 
Mi,M 2 < 10~ 4 M 3 . This is phenomenologically very interesting because it was found that 
under these conditions the effective Majorana mass in neutrinoless double beta decay might 
be related to the CP violating phase controlling leptogenesis Furthermore, since the 
heaviest right-handed neutrino effectively decouple, this scenario might provide a natural 
framework for models with two right-handed neutrinos only 
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The virtue of this approach is that all mass matrices are treated in the same way, pro- 
viding a simple framework which can be easily embedded into GUTs. If embedded into 
simple GUTs, the third generation Yukawa coupling unification (at least approximate) is 
inevitable which is theoretically very appealing and it has interesting consequences for phe- 
nomenology |20|]. Note that in this framework it is achieved without making one generation 
different from others at a fundamental level. On the other hand, the spectrum of the first 
two generations of quarks and charged leptons crucially depends on small perturbations. 
In the neutrino sector, the heavier two neutrinos are model dependent, while the mass of 
the lightest neutrino in this approach does not depend on perturbations in the leading or- 
der. The right-handed neutrino mass scale can be identified with the GUT scale, in which 
case the mass of the lightest neutrino is given as (m^ op /M GUT ) sin 2 6*23 sin 2 6 12 in the limit 
sin 9 13 ~ 0. 

We do not provide any understanding of the origin of universal mass matrices and their 
perturbations. It is not straightforward to construct such models with family symmetries. 
Nevertheless, it is worthwhile to look for alternatives: extra dimensions or composite models. 
No matter what the origin is, having all three generations indistinguishable in leading order 
is certainly something one would like to see in the fundamental theory. After all, the three 
generations have exactly the same quantum numbers in the standard model and even in 
simple GUT models. Why should their Yukawa couplings be so different? 
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